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■ We evaluate the scattering functions of a gas of spin-polarized, non-interacting fermions confined 

(N ■ in a quasi-onedimensional harmonic trap at zero temperature. The main focus is on the inelastic 

, scattering spectrum and on the angular distribution of scattered light from a mesoscopic atomic 

^ ■ cloud as probes of its discrete quantum levels and of its shell structure in this restricted geome- 

pL^ ' try. The dynamic structure factor is calculated and compared with the results of a local-density 

approximation exploiting the spectrum of a one-dimensional homogeneous Fermi gas. The elastic 
and inelastic contributions to the static structure factor are separately evaluated: the inelastic term 
becomes dominant as the momentum transfer increases, masking a small peak at twice the Fermi 
momentum in the elastic term but still reflecting the discrete quantum levels of the cloud. A fast and 
' accurate numerical method is developed to evaluate the full pair distribution function for mesoscopic 

^ , inhomogeneous clouds. 
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I. INTRODUCTION 



There currently is considerable experimental activity in trapping and cooling dilute gases of the ferniionic alkali 
atoms ^"^K [Q and ^Li in the quantum degeneracy regime. A quasi-onedimensional (quasi-lD) and almost ideal 
^ \ neutral Fermi gas can be realized in these experiments by using a very anisotropic, cigar-shaped magnetic trap. In 
Q ■ this arrangement the fermionic atoms can be placed into a single Zeeman sublevel, yielding a spin-polarized state in 
O ' which their mutual interactions are extremely weak [^[^, and are subject to approximately harmonic confinement in 
the axial direction and to tight confinement in the other directions. 
I ' The possibility of experimentally studying such a close realization of an ideal ID Fermi gas in harmonic external 
^ I potential presents great interest for the theory of quantum gases. To begin with, the ID model of non-interacting 
' fermions has the same particle density profile and thermodynamic and dynamic properties as the ID model of point- 
[ like bosons with hard-core interactions [g . Given this wider relevance of the model, harmonic confinement is probably 
■ the simplest way to realize an inhomogeneous quantum fluid on which to test the ideas and the implementations of 
Density Functional Theory Of specific interest in the present context is the kinetic energy functional, which 
can be constructed exactly from the particle density in the case of TV fermions moving independently in a harmonic 
oscillator potential 0. Finally, the shell structure that was noticed in the particle density of the Fermi gas in 3D 
"j'-y is greatly enhanced in lower dimensionality . Very eflScient numerical methods have been developed Ppd]] to 

' evaluate ground-state properties of mesoscopic ID clouds containing large numbers of fermions. 

The present work completes our earlier study of the ideal ID Fermi gas under harmonic confinement 11 1 by 



I ■ evaluating the spectral properties of low-frequency radiation scattered from its quasi-lD realization. Previous studies 
' O ' of the scattering spectrum of the homogeneous ideal Fermi gas [^[^ and of the off-resonant light scattering properties 
^ . of an ideal Fermi gas under isotropic harmonic confinement ]H have referred to the 3D case. Our objective is to 
examine how the shell structure of the particle density profile and the single occupancy of discrete quantum levels 
in this specific quasi-lD geometry are reflected in the light-scattering properties and, at a deeper level, how the shell 
structure appears in the pair distribution function as a gauge of the spatial correlations induced by the Pauli exclusion 
principle between pairs of fermions. 



5_] ■ The layout of the paper is briefly as follows. In Sections ^ and III we recall from the work of Javanainen and 
I Ruostekoski [|l^ the basic facts about light scattering from an atomic cloud and explicitly set out how the inelas- 
tic scattering spectrum and the angular distribution of the scattered light are related to the time-dependent pair 



distribution function. In Sec. IV we present an evaluation of the dynamic structure factor of a mesoscopic Fermi 
gas in harmonic confinement and compare its results with those obtained from a local-density approximation (LDA) 
exploiting the spectrum of the homogeneous gas. Section |^ evaluates both the elastic and the inelastic contribution 



to the static structure factor, again comparing the latter results with the LDA. In Sec. VI we present numerical results 
for the equal-times pair distribution function, which are obtained by an e xtension of the Green's function method 
previously developed for the one-body density matrix. Finally, Sec. VIl summarizes our main conclusions. 
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II. LIGHT SCATTERING FROM AN ATOMIC CLOUD 



We consider a confined cloud of spin-polarized, magnetically trapped atoms in a quasi- ID configuration. We focus 
on a light scattering experiment with the incident beam propagating along the z axis, orthogonal to the x axis which 
is the axial direction of the magnetic trap (see Fig. |l|) . The positive- frequency component of the incident electric field 
is 

E^(^,t) = ^feV('=^^--"*^ (1) 

where £ is the field amplitude, e is the polarization, kh and lol are the light wave number and frequency. In regard to 
the internal states we model each atom as a two-level system, i.e. an excited state |e) and a ground state \g) separated 
by an energy hu>A- This corresponds to a suitable choice of the Zeeman sublevel for the atoms in the ground state 
and, at the same time, of the light polarization. Angular momentum degeneracy of the internal states in the case of 
an unpolarized cloud has been treated by Javanainen and Ruostekoski p^ . 

In the Born approximation the scattered light at a point on the screen designated by the vector R (see again Fig. |l|) 
is given by Q 

E+(R,i)=£:-^^[nx(nxd-)](e.d+) / d\'e-^'^-'i,\v'M{^,t). (2) 
STreofto ti Jcioud 

Here 8 = ljl — oja is the detuning between the laser frequency and the atomic transition, fi = R/i?, d+ = 'Dt\e){g\ 
and d~ = T) €\g){e\ are the raising and lowering parts of the dipole-moment operator of amplitude I?, k = fcifi — 
is the wave vector transfer and r' indicates an atom in the cloud (see Fig. |l|). Finally, V'(r,f) is the annihilation field 
operator for atoms in the internal state \g). 

The spectrum of scattered radiation from atoms in the ground state at the point R is 

X(k,c.) = l-M{e,n)I{R)^ j dt j d^n J e'('^*-''-("-^-^^» (^t (j.^^ t)^(ri, i)V-t(r2, 0)^(r2, 0)) . (3) 

Here M(e, n) = (1 -I- cos^(6'))/2 is the geometric factor for detection at an angle 9 in the case of circularly polarized 
light; I{R) = 'D'^uj\/{2>2iT'^eQC?B?) is the total intensity scattered by a single dipole, and = T>£/{2h) is the Rabi 
frequency. The triple integral on the RHS of Eq. (^) at w 7^ is the dynamic structure factor ^(k, lo) for inelastic 
scattering. 

The angular distribution of the scattered light is 

I(k) = J du;I{k,uj) , (4) 

the angular dependence being implicit in the direction of k. The distribution of scattered light can be decomposed 
into two components: an elastic term, originating from the diffraction of the light by a finite object of given density 
profile, and an inelastic term determined by the excitations of the quantum fluid. 



III. GENERAL DEFINITIONS 

According to Eq. (^ the microscopic quantity which is relevant for the spectrum of scattered light is the time- 
dependent pair distribution function, 

P2(l,2) = (V'^(l)^(l)^t(2)^(2)). (5) 

Here we have set (1) = (ri, ii) and (2) = (r2, t2)- For a non-interacting Fermi gas the expansion of the field operators 
into energy modes Oj, that is '0(1) = ^^p{~'>'£iti)0'i^ yields an explicit expression for the pair distribution 

fmiction in terms of the single-particle orbitals ^^(r): 

P2(l, 2) = ^ 0:(ri)0,(ri)0*(r2)0.(r2)e-^fe-^-)(*-*=)/''/(£.)[l -/(£,)]+ n(l)n(2) 

= S{l,2) + nil)n{2) , (6) 
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where n(l) is the equiUbrium density profile and f{e) is the Fermi distribution function. 
The dynamic structure factor is given by 

Sik,Uj) ^ J dtj d^nj d3^2g^[^t-k.(rl-r2)] 5'(1,2)| 



E 



d3re-^''-</>*(r)^,(r) 



2 



/(£.)[! -/(£,)] 27r<5(c. A). (7) 



By integration over frequency one obtains the inelastic contribution to the static structure factor: 

d'n I d\2 e-^''-('-^-'-^) [P2(ri, rs) - n{v^)n{v2)] . (8) 



Here /92(ri,r2) is the equal-times pair distribution function, which for non- interacting fermions can be expressed as 

P2(ri,r2) = n(ri)(5(ri - - pi(ri, r2)/9i(r2, ri) + n(ri)n(r2) , (9) 

where pi(ri,r2) = {tp'' {ri)tp{r 2)) is the one-body Dirac density matrix. 

The next sections will be devoted to the evaluation of these quantities in a quasi-lD Fermi gas under harmonic 
confinement Vextiic) — mojf^J^x^ -\- A^(j/^ -I- 2^)]/2, that is, a Fermi gas in a cigar-shaped harmonic trap where A 1 
and only the values Uy — Q and = of the transverse quantum numbers are allowed. For simplicity we choose 
h = I, m — 1 and cuho = 1 • 

IV. DYNAMIC STRUCTURE FACTOR 

The transverse excited states are not involved in the excitation processes when (i) the transverse component k± 
of the momentum transfer vanishes, due to orthogonality of harmonic oscillator wave functions, or (ii) the energy 
transfer lo is smaller than the gap between the chemical potential and the first transverse excited state. In these limits 
Eq. reduces to a one-dimensional problem and the transverse-state wave function factorizes out. Following De 
Marco and Jin [|| we evaluate the overlap integral / = J dxi exp{ikxXi)(j)* {xi)(t)j{xi) in momentum space to obtain 

1 = J ^4>*iq + kj2)4>j{q~kx/2), (10) 

where ^i{p) — "Hiip) exp(— p^/2)/7r-'-/^V2*z!. By exploiting the properties of the Hermite polynomials Hiip) [|^, we 
finally obtain the expression for the dynamic structure factor as 

^-1 -I /t.2\'* 

i=max{Af-ft,0} ^ ^ \ / 

Here h is an integer corresponding to a single-atom excitation of h quanta of the harmonic oscillator and L^{x) is 
the i*^ generalized Laguerre polynomial of parameter h. The spectrum in Eq. ( pT| ) is discrete due to the external 
harmonic confinement, but can of course be approximated as a continuum in the limit uj ^ oJho- Figure ^ illustrates 
the spectrum at various numbers of particles, for two values of the transferred wave vector in the x direction. 

The main qualitative features of the spectrum can be understood within a local-density description, which exploits 
the results for the dynamic structure factor Shomikx,ui) of the homogeneous ID system in a medium with a locally 
varying density profile n{x) to predict the spectrum of a strictly ID system as 

SLDA{kx,uj)^ J dxn{x)Shoni{kx,i^; ^J-{x)) ■ (12) 

The spatial dependence of the chemical potential ^ = kp/2 = 7r^n^/2 is determined by the relation fi{x) = ^{n{x)). 
In the case of external harmonic confinement we use the expression for the homogeneous ID Fermi gas, 

if \uj2{kx)\ < UJ < uji{kx) 

Shoni{kx,Uj) = { (13) 

otherwise 
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where 1^1.2 = ^x/^ i k^kp, together with the density profile in LDA, 

n{x) = -\/2N-x^ (14) 

TT 

to obtain the dynamic structure factor: 



SLDA{kx,i^) = -r~ 
kx 

The LDA curves for the dynamic structure factor are also shown in Fig. y, together with the spectrum of the 
homogeneous gas. Evidently, the effect of the confinement is to change very substantially the spectrum with respect 
to the homogeneous case, and the general spectral shape is well described within a local-density picture. 




V. ANGULAR DISTRIBUTION OF THE SCATTERED LIGHT 



The angular distribution of the scattered radiation is related through Eqs. (^) and (Q) to the equal-times pair 
distribution function. Both an elastic term and an inelastic term contribute to its expression: the former is related 
to the particle-density profile of the cloud while the latter arises from the internal excitations of the cloud. 



A. Elastic contribution 



The elastic contribution to the diffraction pattern of a quasi-lD gas is given by 

2 



5e(k) =e 



dxe *''="^ri(a;) 



(16) 



In evaluating this term, which is of order A^^, we have looked for the effect of the shell structure of the density profile 
on the angular distribution of the scattered light. 

With this aim, we have compared the Fourier Transform (FT) n(kx) of the ID density profile with its LDA 
expression, which can be evaluated in terms of the Bessel function Ji as |lO| 

flLDAikx) = -, JliVmkx) . (17) 



kx 

The same behaviour is found for these two curves in the central peak of the diffraction pattern (Fig.^(a)), where the 
dominant contribution is due to the finite size of the cloud. The FT of the exact profile differs instead from that of 
the LDA (Thomas- Fermi) profile in that it contains a side peak at a value kx — I-k j of the transferred momentum, 
where Ax is the average distance between two nearest maxima in the shell structure. An estimate of Ax can be made 
by taking the maxima to be uniformly distributed along the LDA profile, yielding 

2Rlda 2^/2 aho ,^„s 
Ax ~ = = — , (18) 



where Rlda is the axial Thomas Fermi radius of the cloud. For A^=100, 500 and 1000 we have recognized this side 
peak in the diffraction pattern (see Fig.^(b) for A^=100). As noticed earlier jl^] its position corresponds to kx — 2kp. 
However, the relative intensity of the peak is extremely weak and as we expected it decreases with increasing N ^ 
as the density profile gets closer to the Thomas-Fermi one. At higher values of the momentum transfer the elastic 
contribution to the diffraction pattern approximately vanishes, while the Thomas-Fermi pattern has an oscillating tail 
which still reflects the finite extension of the profile in space. 

The total distribution of the scattered light is the sum of the elastic and inelastic contributions: while for small 
values of kx the elastic contribution is dominant but does not distinguish between the true and the Thomas-Fermi 
profile, the inelastic contribution overcomes it at finite values of kx- This is shown immediately below. It is thus very 
difficult to resolve the specific effect of the shell structure by this method. 
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B. Inelastic contribution: the static structure factor 



The inelastic contribution to the intensity of scattered radiation is obtained from the static structure factor in 
Eq. (^). This contribution is proportional to the number N of particles in the sample. In the case of quasi-lD 
harmonic confinement, by the same procedure which has led to Eq. ( pT[ ) we obtain an analytic expression for S'i(k) 
in terms of the generalized Laguerre polynomials: 



J! 



Lrikl/2) 



(19) 



In Figure ^ we show the function 5i(k) at fc^ = for various values of the particle number N. The condition fc^ ~ 
applies to small angles in Fig. |l]. The quantum effect of single-level occupancy in the external potential is most visible 
in the first derivative of Si{kx), as is shown in the inset in Fig. ^. 

The local-density result, obtained by integration of the ID Eq. (|l5|) yields the following expression for a strictly ID 
system: 



^LDA 



2NkrcT&n{kx/^Jm - k^) + k^y/m - k^/A /n if k^ < 2^/2N 



N 




(20) 



It follows from the discussion in Sec. |^that the ID LDA provides a good description of the structure factor at fcj^ = 
as N increases. In fact, although the LDA does not take shell effects into account, it gives the general behavior of 
Si{kx) already for small numbers of fermions {N > 20). The LDA curve has not been shown in Fig. ^ since it is 
indistinguishable from the N = 100 result. 



VI. EQUAL-TIMES PAIR DISTRIBUTION FUNCTION 

Because of its definition as a double space integral in Eq. (H), the inelastic term in the static structure factor displays 
only part of the information which is contained in the equal-times pair distribution function /92(ri,r2) a-s defined in 
Eq. (^). In this section we give a full evaluation of this function for the quasi-lD gas of present interest. Our main 
objective will be to explicitly show how /92(ri,r2) reflects the Pauli exclusion principle and the shell structure of the 
fermion cloud. The same results apply to a ID Bose gas with hard-core interactions, since the boson-fcrmion mapping 
holds for the equal-times pair distribution function |15| . 

With the notation (j)_L{y, z) for the transverse ground-state wave function of the quasi-lD system of non-interacting 
fermions in the external confining potential, we can write the pair distribution function as 

P2(ri,r2) = n(ri)(5(ri - V2) - w{yi^,Y2i_) [F{xi,X2) - n{xi)n{x2)] (21) 

where rj_ = (y, z) and ■w{vi^,V2i_) = |0±(2/i, 0i)P|0_l(2/2, 22)^- The function F{xi,X2) is given by 

N 

F{xi,X2) = 4>*{x^)'^o{x^Wj{x2)Ux2) (22) 

and will be calculated below by an efficient numerical method based on an extension of the Green's function method 
used previously to evaluate particle and kinetic energy densities. Direct calculation can be used for limited 

numbers of fermions [|l5|, but the present method is easily applied to handle mesoscopic clouds containing a large 
number of fermions. 



A. Method 

We first rewrite Eq. ( [2^ ) in terms of the Green's function G{x) — [x — x + ie)^^ in coordinate space, where x is the 
position operator. We obtain 
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1 ^ „ 

F{xi,X2)^^ lim Im(0j|G(xi)|0j) Im(0j|G(x2)|0i 



— ■ lim Tr 



|lm |^GAr(a;i; 



Im 



G 



N 



(23) 



where are the eigenstates of the ID Hamiltonian in the coordinate representation and Gn{x) is the first N x N 
block of the matrix G(x). We then use the property Im^ • Im_B ~ (l/2)[Re(A • B*) — Re{A ■ B)] to obtain 



Fixi,X2) = lim (RcTr [Gjv(a;i) • G^(x2)l - ReTr [GAr(a;i) • Gn{x2)\ } 



(24) 



We can now use the relation TrQ ^ = [d\ogdet{Q + XI)/dX\x=o [Q, with I being the identity matrix N x N, to 
obtain the final expression 



F{xi,X2) 



lim Re-^ 

— hm Re— 

27r"^ e^o+ oX 



logdet (^{x2 ~ C{x2) - ie) ■ {xi - i{xi) + ie) + X. 
logdet (^{x2 - i{x2) + ie) ■ {xi - i{xi) + ie) + X 



x=o 



x=o 



(25) 



Here, ^{x) is a renormalized position operator jl^], reduced to the N x N subspace. 

While Eq. (^5|) holds for any ID confining potential, this approach is particularly useful for the case of harmonic 
confinement where the representation of the position operator is tridiagonal. In this case the renormalized operator 



^(x) reads [^(a 



J2J 



[x]ij if (i,j) (iV,iV) and [^{x)]n,n = xn,n{x), with 



XN,Nix) = 



N/2 



x + le ■ 



(iV+l)/2 ■ 



(26) 



X + le ■ 



The evaluation of F{xi, X2) is thus reduced to the calculation of the determinant of pentadiagonal N x N matrices. 
If we make a partition into blocks Ai and Bi j of dimension 2 x 2, a pentadiagonal matrix Q takes a tridiagonal form. 
For even N we can write (see also ) 



/ Al Bi,2 

i?2,l A2 ^2,3 



Q = 



\ 



B' 



3,2 



V 



(27) 



■ • An/2 I 

and the determinant of such a matrix can be factorizcs into the product of 2 x 2 matrices as 

det Q = ]J det Ai . 



(28) 



Here, A\ = A\ and Ai — Ai — Bi,i^i{Ai^i) ^Bi^i^i for i > 1. If is odd we obtain an analogous expression in which 
the last block of the partition is a 1 x 1 matrix. 



B. Numerical results 



By the method outlined above we have calculated the longitudinal contribution to the equal-times pair distribution 
function, {xi, X2) — n(xi)n(x2) — F{xi,X2), for up to = 100 fermions without particular numerical efforts. Since 
this function presents in general a number of maxima of order 7V^, for the sake of clarity we have shown in Fig. |^ the 
full result only for the case of = 4 fermions. 

The numerical results for higher values of the particle number N are shown in Fig. ^ through a section of the pair 
distribution function taken at the value R = {xi +X2)/2 = of the center-of-mass coordinate. In the confined gas the 

effect of Pauli exclusion is also observable in real space as a depression of P2\xi, X2) at short distance r = {xi — X2)- 
This is more clearly illustrated in the inset of Fig. ^ through an enlargement of the region near r ~ 0. At larger 
distances the function F{xi, X2) goes rapidly to zero and the pair function decouples into the product of two particle 
density profiles. For other sections of the pair function a similar behavior is found with a reduced number of peaks in 
the profile. This is shown in Fig. |^ for the illustrative case of A^ = 4 fermions. 
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VII. SUMMARY AND CONCLUSIONS 



In summary, we have evaluated the scattering functions of quasi-lD mesoscopic clouds of non-interacting fermions 
in an external harmonic potential. A number of conclusions can be drawn from our results. 



The inelastic scattering spectrum of the inhomogeneous gas, displayed in Fig. S, reflects its restricted geometry 
through (i) the form factor coming from the tight transverse confinement, and (n) its general shape, which is re- 
markably different from that of the homogeneous gas, as can be predicted already in a local-density calculation (see 
Eqs. (p^-(p^). In principle the spectral intensity has a discrete structure owing to the level spacings induced by the 
axial confinement, but this spectral structure is rapidly smoothed away as the number N of fermions in the cloud 
increases. 

The elastic contribution to the scattering cross section is dominated by the main peak centered at zero momentum 
transfer and broadened as a consequence of the finite axial size of the cloud (see Fig. || (a)). The shell structure of 
the cloud appears in the elastic scattering intensity only as an extremely weak peak located at momentum transfer 
close to twice the Fermi momentum. This peak lies in the far tail of the central elastic peak (see Fig. I (b)) and will 
be masked by the inelastic scattering contribution to the diffraction pattern. 

The inelastic contribution to the static structure factor in Fig. |^ is obtained by integration of the dynamic structure 
factor over energy transfers. Similar comments apply, therefore, to these two scattering functions: thus, Si{kx) is 
clearly reflecting the shell structure of the cloud only for limited values of the number of fermions (see the inset of 
Fig. U). The general shape of Si{kx) in the main body of Fig. |4|is again related through the local-density scheme to 
the static structure factor of the ID Fermi gas. 

The equal-times pair distribution function contains a much greater and detailed amount of information on the 
relative spatial distribution of pairs of particles in the inhomogeneous Fermi gas. Here it is, of course, a function 
both of the relative distance r of the pair and of the coordinate R of its center-of-mass along the x axis. The Pauli 
"exchange" hole surrounding each fermion is a well-known feature of the pair distribution function in homogeneous 
Fermi gases and is a very prominent feature in the spin-polarized inhomogeneous system. It is seen in Fig. |^ as the 
deep valley running in a direction parallel to the R axis and cutting into two separate parts the contour plot at r = 0. 
It is again seen in the inset of Fig. ^ at i? = and in Fig. ^ at various values of R. In addition to this basic feature, 
the two-body distribution displays a great deal of secondary structures, of order as already noted, which descend 
from the combination of the Pauli principle and of the axial confinement. 

In conclusion, the single occupancy of the discrete quantum levels in the quasi-lD Fermi gas may become accessible 
in light scattering experiments by looking at the inelastic scattering spectrum. It also affects to some extent the 
angular distribution of the inelastically scattered light. The elastic contribution to the scattering is quite negligible 
in this respect. 



A.M. acknowledges useful discussions with Dr. I. Carusotto. This work was partially supported by MURST under 
the PRIN2000 Initiative. 
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0.05 0.1 0.15 0.2 0.25 2 4 6 8 

FIG. 2. Dynamic structure factor S{kx,uj) of a ID harmonically confined Fermi gas (in units of (jjJ^^) as a function of huj/Ep, 
where Ef = NhuJho- Left panel: at kx = 0.1 kp and numbers of fermions A'^=500 (solid line) and A'^=1000 (solid bold line). 
Right panel: at fci = 2 kp and N=2Q (solid line) and iV=40 (solid bold line). In both panels the LDA spectrum (dotted 
line) and the dynamic structure factor of a homogeneous ID Fermi gas (dashed hnc) arc plotted in the same units. The exact 
S{kx,ijo) is shown for clarity as constructed from discrete frequencies by joining vertical lines having height proportional to the 
oscillator strength. 



1.5x10 




V2kp(N) ky2k,(N) 
FIG. 3. Central peak (left panel) and side peak (right panel) of the elastic contribution to the static structure factor, 

normalized to N = 100 fermions, as a function of the wave number kx (in units of 2kF{N) = 2^/2N aj^^) at fixed k± = 0. Both 
the true profile (continuous line) and the LDA profile (dotted line) are shown. 
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ky2k,(N) 

FIG. 4. Static structure factor for a quasi-lD Fermi gas in harmonic confinement as a function of the wave number kx (in 
units of 2kFiN)) at k± = 0. The inset shows the first derivative of Si{kx) in the same reduced units. The numbers of fermions 
are N =4 (short-dashed line), N =20 (long-dashed line) and N =100 (solid line). 




FIG. 5. Longitudinal contribution p2\r; R) to the pair distribution function for N=4: fermions in quasi-lD harmonic con- 
finement, in units of a^^, as a function of the center-of-mass coordinate R and of the relative coordinate r, both in units of 
a/to- 
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10 20 30 

FIG. 6. Section of p2\r\R) at _R = for A'^=20, 50 and 100 fermions at _R = (from bottom to top), in units of a^^, as 
a function of the relative coordinate r/a^o of the pair. The inset shows an enlargement of the region near r = 0, in the same 
units. 




2 4 6 

r/aho 



FIG. 7. Sections of p2 {r;R) for values of R—Q (solid line), R = a^o (dashed line) and R — 2aho (dotted line) at N=A. 
Notations and units are as in Fig. M. 
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